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Abstract. A complete set of d + 1 mutually unbiased bases exists in a Hilbert 
spaces of dimension d, whenever d is a power of a prime. We discuss a simple 
construction of d+ 1 disjoint classes (each one having d— 1 commuting operators) 
such that the corresponding eigenstates form sets of unbiased bases. Such a 
construction works properly for prime dimension. We investigate an alternative 
construction in which the real numbers that label the classes are replaced by a 
finite field having d elements. One of these classes is diagonal, and can be mapped 
to cyclic operators by means of the finite Fourier transform, which allows one to 
understand complementarity in a similar way as for the position-momentum pair 
in standard quantum mechanics. The relevant examples of two and three qubits 
and two qutrits are discussed in detail. 



PACS numbers: 03.65.Ta, 03.65.Ca, 03.65.Ud, 03.67.-a 
1. Introduction 

The concept of complementarity is a direct nontrivial consequence of the superposition 
principle and distinguishes purely quantum systems from those that may be accurately 
treated classically. Therefore, a thorough understanding of this idea is of fundamental 
importance for a correct interpretation of quantum mechanics [1]. 

In short, Bohr idea of complementarity could be loosely formulated by stating 
that, in order to understand a quantum phenomenon completely we need a 
combination of mutually exclusive properties: the precise knowledge of one of them 
implies that all possible outcomes in the other are equally probable. 

Perhaps the best textbook illustration of complementarity is that the observation 
of interference and the knowledge of the path followed by the interfering particle are 
mutually exclusive. This is often expressed as the statement that all quantum objects 
exhibit particle-like or wave-like behavior under different experimental conditions. 
Standard examples, such as Einstein recoiling slit [2], Feynman light-scattering 
arrangement [3] or Heisenberg microscope [4] , are always explained in terms of position 
and momentum. It is thus natural that links between canonical conjugacy and 
complementarity have been fully explored [5-7]. 

However, most if not all of the recent examples of complementarity involve 
finite-dimensional systems (for a complete and up-to-date review see [8]). Note 
that, contrary to what one would expect, the finite dimensionality of the system 
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space introduces difficulties when dealing with uncertainty relations, because the 
commutators between complementary observables are operators instead of c- numbers, 
which makes the analysis of the problem more involved. The situation has been 
recently addressed for two-dimensional spaces [9]; however, there are subtle aspects 
that cannot be fully encompassed by these easy-to-understand systems. 

In finite-dimensional systems, complementarity is tantamount to unbiased- 
ness [10,11]: each eigenstate of any measurement is an equal- magnitude superposition 
of the eigenstates of any of the complementary measurements. This leads naturally 
to the concept of mutually unbiased bases (MUBs), which have recently been con- 
sidered with an increasing interest because of the central role they play not only 
in understanding complementarity [12-14], but also in specific quantum information 
tasks, such as protocols of quantum cryptography [15,16], Wigner functions in discrete 
phase spaces [17,18], or the so-called Mean King problem [19-23]. 

For a d-dimcnsional system it has been found that the maximum number of 
MUBs cannot be greater than d + 1 and this limit is reached if d is prime [24] or 
power of prime [25,26]. Remarkably though, there is no known answer for any other 
values of d, not even for d — 6. Recent works have suggested that the answer to 
this question may well be related with the non-existence of finite projective planes of 
certain orders [27, 28] or with the problem of mutually orthogonal Latin squares in 
combinatorics [29,30]. 

Quite recently, a number of papers have addressed the explicit construction 
of MUBs for dimensions that are prime or composite (i. e., power of a prime), 
exploiting different algebraic properties [31-36]. In this paper we give an explicit 
construction with a different method that resorts to elementary notions of finite 
field theory. We wish to emphasize the distinct features of our approach: first, we 
recall that complementarity for the position-momentum pair is implemented by the 
Fourier transform, which exchanges both operators. Therefore, we construct classes 
of maximally commuting operators and map them using the finite Fourier transform 
and an additional diagonal operator. In consequence, we obtain in a systematic way 
the whole family of complementary operators and not merely MUBs. Additionally, 
the final expression for these MUBs is compact and can be immediately expressed in 
different bases, in some of which they appear as tensor products of generalized Pauli 
matrices. In summary, we hope that our unified construction provides a simple picture 
of complementarity for both prime and composite dimensions. 

2. Complementary operators in prime dimension 

We consider a system living in a Hilbert space Tid, whose dimension d is a prime 
number. It is useful to choose a computational basis \n) (where n = 0, . . . , d — 1) in 
Tid and introduce the basic operators 

X\n) = \n + 1), 

(2.1) 

Z\n) =w n \n), 

where 

to = exp(2m/d) (2.2) 

is a dth root of the unity and addition and multiplication must be understood modulo 
d. These operators X and Z, which are generalizations of the Pauli matrices, were 
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studied by Patera and Zassenhaus [37] in connection with additive quantum numbers, 
and have been used recently by many authors in a variety of applications [38-40]. 
They generate a group under multiplication known as the generalized Pauli group and 
obey 

ZX = luXZ, (2.3) 

which is the finite-dimensional version of the Weyl form of the commutation relations. 

According to the ideas in [31], we can find d+ 1 disjoint classes (each one having 
d— 1 commuting operators) such that the corresponding eigenstates form sets of MUBs. 
The explicit construction starts with the following sets of operators: 

{Z k }, fc = l,...,d-l, 



(2.4) 



{{XZ m ) k }, fc = l,...,d-l, m = 0,...,d-l. 
One can easily check that 

Tr(Z fe Z fe 't) = dS kk >, Tr(X k X k '^) = dS kk >, 



(2.5) 



Tr[{XZ m ) k {XZ m ') k '^] = d5 kk ,6 mm ,. 

These pairwise orthogonality relations indicate that, for every value of m, we generate 
a maximal set of d — 1 commuting operators and that all these classes are disjoint. 
In addition, the common eigenstates of each class m form different sets of unbiased 
bases. We shall refer to these classes as multicomplementary. 

We would now like to make the very important observation that, starting from 
Z, it is possible to obtain any element of the form (XZ m ) k by using a combination of 
only two operators F and V defined as follows: F is the finite Fourier transform [41] 

1 d-1 

F=- / = £ u nn '\n)(n>\, (2.6) 

n.n'— 

and V is the diagonal transformation (assuming d is odd) 

d-1 

V = J2u- (n2 - n){d+1)/2 \n)(n\. (2.7) 

n=0 

Indeed this is the case, since one easily verifies that 

X = F*ZF, (2.8) 

much in the spirit of the standard way of looking at complementary variables in the 
infinite-dimensional Hilbert space: the position and momentum eigenstates are Fourier 
transform one of the other. On the other hand, the diagonal transformation V acts 
as a Z-right shift: 

XZ rn = V frn XV m . (2.9) 

The case d = 2 needs minor modifications. In fact, it turns out that one cannot 
find a diagonal unitary transformation V such that X — ► XZ. For this reason, instead 
of XY the matrix Y is defined as iXZ, so that Y = X V, where V is 

v-(ll). <"») 

The construction of multicomplementary operators is otherwise identical to the case 
where d is an odd prime. 
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3. Complementary operators in composite dimensions 

3.1. Constructing multicomplementary operators 

For all its simplicity, the construction of the previous Section fails if the dimension of 
the system is a power of a prime. A simple illustration of this is obtained in dimension 
4 = 2 2 . According to equation (|2.1(l . the operators X and Z are simply 



X = 



(° 


i 













1 














1 


V i 








o J 



z = 






-1 





(3.1) 



/ 



Then, for instance, X 2 = —(XZ 2 ) 2 and (XZ 3 ) is proportional to XZ, so operators 
constructed following (|2 .4|> no longer form disjoint sets. The root of this failure can 
be traced to the fact that Z4, the set of integers modulo 4, does not form an algebraic 
field. The same failure generally occurs for any composite dimension d = p n , where 
p is a prime and n is an integer. In short, the construction of multicomplementary 
operators cannot proceed by simply taking powers of some basic elements. 

However, we know there exists (up to isomorphisms) exactly one field, written as 
F<j, with d elements when d = p n . If d = p is prime, the field essentially coincides with 
Z p . We briefly recall the minimum background needed to proceed. For more details 
the reader is referred to the pertinent literature [42]. The field can be represented 
as the field of equivalence classes of polynomials whose coefficients belong to Z p . The 
product in the multiplicative group F5 (i. e, excluding the zero) is defined as the 
product of the corresponding polynomials modulo a primitive polynomial of degree n 
irreducible in 1 p . In fact, ¥* d is a cyclic group of order d— 1: it is generated by powers 
of a primitive element a, which is a monic irreducible polynomial of degree n. This 
establishes a natural order for the field elements, and we use this order to label the 
elements of a basis in TLd as follows: 



-d 



{|0)»,|a%..,|a^)}. (3.2) 

to the problem of MUBs in composite dimension consists in using 
instead of natural numbers, to label the classes of complementary 



Our solution 
elements of 
operators. 

Next we define the trace of a field element 9 £ 



as 



tr(0) = + B v + 6 P ~ + . . . + 9 P . (3.3) 

Note that we distinguish it from the trace of an operator by the lower case tr. The 
trace has remarkably simple properties, the most important for us being that it is 
linear and that it is always an element of the prime field Z p . 

For the additive group in the field F^ we can introduce additive characters as a 
map that fulfills 

x(0t)x(02)=x(0i + 6 2 ), 0i,02 eF rf . (3.4) 

All of these additive characters have the form 

'2m 

(3.5) 



X(0) = exp — tr(0) 
L P 

For future reference, we quote the property 

E *w = °. 



(3.6) 
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which leads to the relation 

d-2 

=d6 e ,o-l. (3.7) 

k=0 

We start by introducing diagonal operators with respect to the basis (|3.2fl as 
follows: 

d-l 

Z q = \0)(0\ + J2x(a q+k )\a k )(a k \, q = 0, . . . ,d - 1. (3.8) 

fc=i 

This definition implies 

Z q \a k )= X (a q+k )\a k ), (3.9) 
and the combination property 

d-l 

Z q Z q , = Z (q)+(ql) = \Q){Q\+Y.x{a q+k +a q ' +k )\a k ){a k \. (3.10) 

fc=i 

These are quite natural generalizations of the properties of matrices in the class {Z k } 
in (12.41) . since the \a k ) are eigenstates of Z q . 

In a similar fashion the operators X q are defined as 
d-l 

X q = ^\a k + a q )(a k \ + \a q )(Q\, q = 0,...,d-2, (3.11) 

fc=i 

so that they act as operators shifting \a k ) to \a k + a q ) and satisfy the combination 
rule 

d-l 

X q X q , = X (q)+(q ,) = \a" + a"')(0\ + ]T|« fc + a" + a"')(a k \. (3.12) 

k=l 

Because elements of the field close under addition, a k + cfl is another element in the 
field: it must be that there is some number L(n), called the Jacobi logarithm, such 
that 

a k + a q = a k+L{q - k \ (3.13) 

whenever a k + ofl ^ 0. In applications, a k + ofl can be found from the summation 
table of ¥ d . 

The finite Fourier transform F, when expressed in terms of the basis \a k ), takes 
the form 



d-l d-l 



(3.14) 



|0)(0|+ x(a k ' +k )\a k '){a k \+J2(\0){a k \ + \a k ){0\) 

k,k' = l k=l 

This definition satisfies the natural property 

F t F = FF t =J. (3.15) 
The operator F is set up to transform the operators Z q into X q : 

X q =F^Z q F. (3.16) 
The Weyl form of the canonical commutation relations is now 

Z q X q , = X (a q+q ')X q ,Z q . (3.17) 
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The operators Z q and X q have been designed to be d- "periodic" , in the sense that 

Za = Zo, Xd — Xq. (3.18) 

In fully analogy with the sets in (|2.4[1 . we can generate operators from X q and Z q ; 
they will be of the form X q Z r . Linear independence and orthogonality are guaranteed, 
in the sense that [compare equation 1)2.5(1 ] 

Ti(Z q z\, ) = d 6 qq , , Tv(X q xl, ) = d Sqq> , 

(3.19) 

Tr{(X q Z r )(X q rZ r 'y] =dS qq '5rr'- 
The commutation relations read 

[X q Z r ,X q ,Z r ,] = X (g)+(s , ) Z (r)+(r , ) [ X (a«'+ r ) - X (a q+r ')}. (3.20) 
It is clear from (|3.19|1 and l|3.2U|l that the sets [compare equation (|2.4p ] 
{Z q }, q = 0,...,d-2, 

(3.21) 

{X q Z q+r }, q,r = 0, . . . ,d- 2, 

are disjoint and that every element of a set with a fixed value r commutes with every 
other element in the same set: they define multicomplementary operators. 

Finally, let us consider the form of the diagonal operators similar to i|2.9[) 
transforming X q to X q Z r . If we restrict to odd dimensions, we have 

d-i 

^ (r) = |0)(0|+^x(2- 1 «'' +2fc - 9 )| a fc )(a fc |, q,r = 0,...,d-2, (3.22) 
k=l 

where x means conjugate character and 2 _1 is an element of Z p ; in particular, if 
p = 2N + 1 we have 2 _1 = N + 1. In this way, one can check that 

vj$JX q vj$. = x^a^W.+r. (3.23) 

Using H3.14|l and (|3.22|) we can generate all the complementary bases. Indeed, 
if the vectors \a) q (q = 0, . . . , d — 2) are the eigenstates of Z , then the whole set of 
complementary bases can be obtained as follows: 

Vf>tFt|S> 9 . (3-24) 
3. 2. Complementary operators for two qubits 

We illustrate our approach with the simplest case of a quantum system of composite 
dimension: two qubits described in a four-dimensional Hilbert space Ti.^. To construct 
multicomplementary operators, we start from the field F4 containing four elements. 
The polynomial 

6 2 +9 + l = (3.25) 

is irreducible in Z2 and the primitive element a is defined as a root of l|3.25|l . In 
consequence the four elements of F4 as in l|3.2|l can be written as 

{0,l,a,a + l}, (3.26) 

where we have taken into account arithmetic modulo 2 and the fact that if a satisfies 
equation (|3.25(l . then we have the relations 

a 2 =a + l, a 3 = 1. (3.27) 
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A direct application of the definition (|3.5[) gives 

X (0) = 1, x (a) = -1, X (« 2 ) = -1, 
Using a the representation where 



|0) = 



X (a 3 ) = 1. 



t 1 ^ 












(°\ 







1 















, \a) = 





, l« 2 ) = 


1 


, \a 3 ) = 





V o J 




V o > 




{0 J 




V 1 / 



(3.28) 



(3.29) 



the matrices Z„ are 



Zo 



/I 



Z 2 = 





-1 






f 









-1 






-1 








1 / 

°\ 




1 / 
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1 








°\ 







-1 
















1 
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-1 / 



The matrix realization of the Fourier transform is 



F = 



and the matrices X a 



X 





/I 
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1 




-1 
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-1 




V 1 




-f 
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°\ 













1 




f 













V 


f 





o J 



1\ 

-1 
-1 

1 J 



1 



V 



XqZq 
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(3.30) 



(3.31) 




1 



(3.32) 



X. 



The rest of the sets are routinely obtained according 13.2111 . 
quote all of them for this example: 



For completeness we 
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-1 
o / 



(3.33) 
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XqZ± — 



X2Z0 — 



XqZ 2 — 



X2Z1 
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(3.34) 



\ 



1 

-1 0/ 



(3.35) 



3. 3. Complementary operators for three qubits 

Our next example is the case of three qubits, the Hilbert space of which is eight 
dimensional. For the field F§ the primitive element a is a root of the following 
irreducible polynomial on F2 

6 3 + 8 + 1 = 0. 

In consequence, the elements of Fg according to l|2.4|l are 



i,or 



1}, 



1. 



= 1. 



{0,1, a, a ,a + l,a +a,a + 
where we have taken into account that 
a 3 = a + 1, a 4 = a 2 + a, a 5 = a 2 + a + 1, a 6 = a 
One can obtain again the additive characters in a straightforward way 
X (0) = l, *(<*) = 1, x(a 2 ) = l, X (a 3 ) = -1, 
X (a 4 ) = l, x(a 5 ) = -l, x(a 6 ) = ~l, x{a 7 ) = -I. 
Using the basis labeling as in (|3.29|l the matrix of the Fourier transform is 



/ 1 
1 



F = 



1 
1 
1 
1 
1 

V 1 

6) are 



1 
1 

-1 
1 

-1 
-1 
-1 
1 



1 
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-1 
-1 
-1 
1 
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-1 
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-1 
-1 
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1 

-1 
1 
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-1 
-1 
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-1 
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-1 
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-1 
1 
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-1 
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-1 
-1 



1 \ 
1 
1 

-1 
1 

-1 

-1 

-1 / 



The matrices Z q (q = 0, . 
Z = diag(l, 1,1, -1,1, -1,-1,-1), 
Z 2 = diag(l, -1,1, -1,-1, -1,1,1), 
Z 4 = diag(l, -1,-1, -1,1, 1,-1,1), 
Z 6 = diag(l, -1, 1, 1, 



(3.36) 
(3.37) 
(3.38) 

(3.39) 



(3.40) 



Z 1 = diag(l ! 1,-1, 1,-1, -1,-1,1), 
Z 3 = diag(l, 1,-1, -1,-1, 1,1,-1), 
Z 5 = diag(l, -1,-1, 1,1, -1,1,-1), 
-1,1,-1,-1), (3.41) 
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and, for instance, 

/0000000 l\ 

1 

1 

1 

1 

1 

1 
\ 1 0000000/ 

The rest of them can be easily worked out. 



Xn = 



(3.42) 



3.4- Complementary operators for two qutrits 

Our final example is the case of two qutrits, described in the nine-dimensional Hilbert 
space Tig. The primitive element of the field Fg is a root of the following irreducible 
polynomial on F3 



9 2 +6 + 2 = 0, 
so that the elements of Fg are 

{0, 1,2, a, 2a + 1,2a- 
and the additive characters are 



2, 2, 2a, a + 2, a + 1,}, 



X(0) = 1, x(a)=w, x (a 2 ) = l, x (a 3 ) = u>, X (« 4 ) = • 
X (a 5 ) = u, x(a 6 ) = 1, X (a 7 ) = w, x(a 8 ) = Q, 



(3.43) 



(3.44) 



(3.45) 



where lu = e 27 ™/ 3 and the bar denotes complex conjugation. Using the same basis as 
(|3.29|l . we have that the matrix of the Fourier transform is 



F=~ 
3 



The operators Z q take the form 

Z Q = diag(l, Q, 1, Q, u>, u, 1, w, Co), 
Z 2 = diag(l, Cb, lu, lu, 1, to, Co, Q, 1), 
Z\ = diag(l, u>, 1, w, w, a), 1, a), w), 
Z 6 = diag(l, w, w, ui, l,Q,w,u}, 1), 
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(3.46) 



Zi = diag(l, 1,uj,lu,lu, l,ui,ui,u>), 
Z 3 = diag(l,w,w, 1,lj,l),u), 
Z§ = diag(l, 1,lu,lu,uj,1,lu,lu,lu), 
Z'j = diag(l, lu,lu,1,lu,lu,lu,1,lu). 



(3.47) 



Note that Z q (g = 1 , . . . , 7) are obtained from Zq by a cyclic permutation of the 
diagonal elements, except the first element that always remains 1. The X q are also 
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easily constructed, and we have, e. g. 
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The diagonal operators V q have the form 
V = diag(l, l,u, l,ui, l,u, I, Co), 



V 1 
V 3 

v 5 

V 7 



(3.49) 



V 2 = diag(l,cj,l,a),l,w,l,a;,l), 
V 4 = diag(l, l,u>, I, Co, l,w), 
Vq = diag(l, Co, 1, lo, 1, Co, 1, lo, 1), 

Note that all the complementary bases can be obtained directly from l|3.24Jl : i. e. 
applying V q F to the basis (|3.44|) . 



diag(l, LO, LO,LO,LO,LO,LO,LO,Lo), 

diag(l, lo, lo, Co, Co, lo, lo, Co, Co), 
diag(l, lo, Co, Co, lo, lo, Co, Co, lo), 
diag(l, Co, Co, lo, lo, Co, Co, lo, lo). 



(3.48) 



4. Multicomplemntary operators as tensor products 

One can establish an isomorphism between the form l|3.21l) of complementary operators 
and its representation in terms of direct product of generalized Pauli operators (|2.1(l . 
This isomorphism can be put forward by showing a one-to-one correspondence between 
the basis H3.2J) and the coefficients of the expansion of the powers of the primitive 
element on, for instance, the polynomial basis, formed by (1, a, a 2 , ... , a™ -1 ). In this 
way we get 



where 



a fc -(cf,c«, 



i), 



4 k) G Z. p 



(4.1) 



E 

1=0 



This allows to rewrite the basis ()3.2ll in the equivalent form 



{\0),\c {k) ,c[ k) , 



{\0),\c { o k) )\c[ k) ) 



(4.2) 



(4.3) 



and the representation of l|3.21|l as a tensor product is now possible. This is due to 
the fact that is isomorphic to Z p x . . . x 7L V , with n products. 

It is worth noting that this isomorphism can be settled also in other bases. For 
example, the so-called normal basis, obtained finding an element (3 G such that 



(4.4) 



is a basis of Fd, is quite useful in applications, since squaring a field element can be 
easily accomplished by a right cyclic shift. Of course, different bases lead to different 

(k) 

coefficients c\ , different sets of commuting operators and different factorizations of 
these operators in tensor products. 
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4-1. Two qubits 

To represent (|3.3U[) and l|3.32[l as a tensor product we follow the described above and 
write the basis (|3.2ti|) in the form (|4.3[l . This yields 

|0> = |00>, |a) = |01), |« 2 > = |11), |a 3 ) = |10). 



(4.5) 



Then, 



Zo 
Zi 

z 2 



(4.6) 



where Z = |0)(0| 



X 
X t 
X 2 



(4.7) 



where X = |0)(1| 



|00)(00|- |01)(01|- |11)(11| + |10)(10| = I®Z, 

|00)(00| - |01)(01| + |11)(11| - 1 10) <10| = Z®Z, 

|00)(00| + |01)(01| - |11)(11| - 1 10) <10| =Z®I, 

|1)(1|. In a similar way the representation of the Xk is: 

|11)(01| + |01)(11| + |00)(10| + 1 10) (00| =X®I, 

|00)(01| + |10)(11| + |11) <10| + |01)(00| = I® x, 

1 10} (01| + |00)(11| + |01)(10| + 1 11) <00| = x®x, 

|1)(0|, so that ZX = -XZ. 
Note that such an asymmetrical correspondence is a result of the Fourier 
transform (|3.31[1 . which is not factorized into a direct product of two Fourier operators. 
The other complementary operators are obtained as a simple product of (|4.6|l and (|4.7() : 

(x Za,x 1 z u x 2 z 2 ) ^{x®z,z®y,y®x), 
{x z 1 ,x 1 z 2 ,x 2 z ) ^ (y®z,z®x 7 x®y), 
(x z 2 ,x 1 z ,x 2 z 1 )^(y®i,i®y,y®y), 

where y — XZ. 

We complete this example using the normal basis {/3,/3 2 }, with a 
f3 2 ,a 3 = f3 + P 2 . Then 

|0) = |00), |a) = |10), |a 2 ) = |01), 

and the new tensor product representatives are 

Zq i — > Z <Si Z, X X <g> X, 

Zn-tZtSl, Xit->X®I, 
Z 2 ^I<g>Z, X 2 i-> I ® x , 

while the rest of complementary operators are given by 

(x z ,x 1 z 1 ,x 2 z 2 )^(y®y,y®i,i®y), 
(x Q z 1 ,x 1 z 2 ,x 2 z Q ) » (y®z,z® x,x®y), 
(x z 2 ,x 1 z ,x 2 z 1 ) ^ {z® x,y® x,x®z). 

The symmetrical aspect of (I4.10f) is a consequence of the factorization of Fourier 
transform in this basis. In fact, we have 



= 111) 



(4.8) 

p,a 2 = 
(4.9) 

(4.10) 
(4.11) 



F = F 2 ®F 2 , 



with 



F 2 



V2 



(|0)(0| + |0)(1| + |1)(0| 



1X1 



(4.12) 



(4.13) 
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4-2. Three qubits 



Mapping to tensor product form is obtained for this case in a manner similar to F4, by 
establishing a correspondence between the bases ll,i.2f and l|4.3|l . Taking into account 
H3.37f) we obtain (by expanding over the polynomial basis {I, a, a 2 }): 



|0) = |000), \a) = |010), 

|a 4 )H011), \<X>) = \\\\), 
In this way one gets 

Z Q h-> Z ® I ® I, 
Zi^>I®I®Z, 
Z 2 ^ I®Z®I, 
Z 3 h^ Z®I®Z, 

z±^ i ® z® z, 

2 5 h Z® Z® Z, 
Z 6 t-+Z®Z®I, 



2 > = |ooi). 
l« 6 ) = |ioi), 



W) 



|110), 
= |100). 



(4.14) 



X 
X 1 

x 2 
x 3 

Xi 

x 5 

Xe 



X 

Id 
Id 
X 
Id 
X 
X 



lid 

Xd 

I® 

$ X 
Xd 
$ X 

$Id 



I 

)I 

X 

9i; 

> X ; 
S> X; 

> X, 



where the operators Z and X are those of the two-qubit example. 
The other sets of commutative operators are as follows: 



{X q Z q } 

jZ q+ 4} 
{X q Z q+5 } 
{X q Z q+6 } 

where q = 



{X q Z q+1 } 
{X q Z q+2 \ 
>- q Z q+3 } 

{XgZ q+ 4) 



(4.15) 



(yn,ixz,izx,yxz,iyy,yyy,yzx), 
(xiz,iyi,ziy,xyz,zyy,yyx,yix), 
(x zi,z x z,i zy,yy z,zy x,y x x,x iy), 
(yiz,iyz,zzy,yyi,zxx,xxy,xzx) : (4.16) 
(xzz,zyz,zzx,yxi,ixy,xyx,yiy), 
(yzz,zyi,zix,xxz,iyx,yxy,xzy), 
(yzi,zxi,ny,xyi,zxy,xyy,yzy), 

, 6 and we have omitted the symbol ® to simplify the writing. Again 
the Fourier transform is not factorized in this basis, but it is possible to find a different 
basis where F factorizes. 



4-3. Two qutrits 

The representation of Z q and X q operators in terms of tensor product is obtained 
using the expansion in the polynomial basis {l,a}: 

|0) = |00), |«> = ]01>, |a 2 )H!2), l« 3 ) = |22), |a 4 ) = |20), 



which leads to 



where now 



* 5 ) = |02), |a 6 ) = |21), |a 7 ) = |ll), |a s ) = |10), (4.17) 

Z ^Z 2 Z 2 ; X ^XI; 
Z } i > Z 2 I, X, i — > I X\ 



Z 2 ^IZ 2 , x 2 ^xx 2 - 

Z Z ^Z 2 Z, X 3 ^X 2 X 2 , 

Z 4 ^ZZ, X 4 ^X 2 I, 

Z$^ZI, X b ^IX 2 ; 

Zq^IZ, X 6 ^X 2 X, 

z 7 ^zz 2 , x 7 ^xx, 



(4.18) 



Z = |0>(0|+w|l)(l|+O|2)(2|, 
A-=|l)(0| + |2)(l| + |0)(2|, 



(4.19) 
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so (|2.1f) holds, since ZX = uiXZ. 

The sets of commutative operators are 

{x q z q } h+ {wz 2 ,z 2 x,xy 2 ,y 2 w 2 ,w 2 z,zx 2 ,x 2 y,yw), 

{X q Z q+1 } i > (W/,/W,WW 2 ,W 2 W 2 ,W 2 /,/W 2 ,W 2 >V,W>V), 
{X g Z g+2 } ^ z 2 , z 2 j; J W 2 , W 2 A' 2 , ^ 2 Z, Z ^ 2 J 2 w, w 

{x 9 z g+3 } ~ (wz,z];jA' 2 ^ 2 }v 2 ,w 2 z 2 I z 2 y 2 j 2 A' 1 iw), . 

{x q z q+A \ ^ (yz, Z ^, Z W 2 , W 2 3^ 2 J 2 Z 2 , Z 2 ^ 2 ,^ 2 W, W)'), 1 Uj 

{x q z q+5 } ^ (yi,iy,yy 2 ,y 2 y 2 ,y 2 i,iy 2 ,y 2 y,yy), 
{x q z q+6 } ^ (x z,zw,wy 2 ,y 2 x 2 ,x 2 z 2 ,z 2 w 2 ,w 2 y,y x), 
{x q z q+7 } ^ {yz 2 ,z 2 w 1 wx 2 1 x 2 y 2 ,y 2 z,zw 2 1 w 2 x), 

where q = 0, . . . , 7, y = XZ, W = A"Z 2 and we have omitted global phases that 
appear in the product of operators. Note that the above sets of commuting operators 
are different from those listed in [34] . 



5. Concluding remarks 

In this paper we have solved the problem of existence and construction of sets of 
MUBs in composite dimensions. Inspired by the the approach developed in [31], 
in which an interesting explicit construction was shown for prime dimension, we 
have used algebraic field extensions to produce a solution for composite dimensions. 
Although other constructive algorithms for solving the MUBs problem in composite 
dimension have appeared, the one presented here does not resort to dual basis, and 
so is completely analogous to the prime-dimensional case. We have also provided a 
simple scheme to cast the sets of MUBs observables as tensor products of generalized 
Pauli matrices. This tensor products can be quite different when we represent the 
finite field in different bases. 

Another major advantage of our approach relies on the use of the finite 
Fourier transform and a diagonal shift operator as maps between maximal classes 
of commuting operators. This is in full agreement with our understanding of 
complementarity in the infinite-dimensional case. 

Of course, there are still open questions: the non-prime dimensional case is a 
challenging issue, or what happens in the limit of high dimensions. In any case, 
the picture developed in this paper is a valuable tool to deal with concepts such as 
entanglement or separability in finite-dimensional systems. 
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